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We employ the previously published neutron scattering experiments[lH3| to suggest that Er2Ti2 07 
has a broken-symmetry multiple-q state with tetragonal magnetic symmetry The ordered mo- 
ments do not appear to lie close to the crystal-field anticipated directions[3] and we suggest that 
the low energy gapless mode, visible in specific- heat measurements [^, is not of the usual transverse 
isotropic Goldstone-mode type, but is longitudinal and is associated with the internal transfer of 
magnetism between distinct magnetic Bragg spots. 



PACS numbers: 03.75.Lm, D75.25.+Z, 75.30.Ds 



I. INTRODUCTION 

There is currently a lot of experimental and theoret- 
ical interest in rare-earth pyrochlore magnets of stoi- 
chiometry, R2M2O7, where R is a rare-earth atom and 
M is an atom that has an available M^"*" state, such as 
tin or titanium. Much of the interest has focused on 
'spin-ice'0, the Ho2Ti207 and Dy2Ti207 compounds, 
due to the proposal that they can be considered to be 
analogue 'monopole' systems [8|. The original motivation 
for studying these systems was the macroscopic degener- 
acy expected for the isotropic Heisenber g ni odcl|9l. The 
compounds Gd2Ti207[l3| and Gd2Sn207Tl| are systems 
that are expected to conform to this model at inter- 
mediate temperature. The issue of how this residual 
degeneracy is lifted at very low temperature then sur- 
faced and further work understanding the dipolar inter- 
action was then relevant [l^. Spin-ice is a system with a 
strong crystal-field interaction forcing the spins to point 
along the natural crystallographic directions whereas the 
gadolinium compounds have essentially no anticipated 
crystal field (subject to some unexpected ESR issues [isll). 
To complete the set it would be nice to have a material 
where the local crystal-field promoted spins perpendic- 
ular to the natural crystallographic directions, and the 
compound Er2Ti207 was thought to play this role[T^. 
Theoretical crystal- field calculations [sj, neutron scatter- 
ing experiments [ll-Q and specific heat measurements j^] 
were all brought forward to promote this picture. Indeed, 
a subtle 'order from disorder' calculationll5| appeared to 
predict the observed experiments. In this article we of- 
fer an alternative view that the spins are oriented quite 
close to the natural crystallographic axes, but are actu- 
ally slightly distorted away into a point-symmetry broken 
orientation. This loss of point-symmetry then leads to a 
continuous degeneracy which explains the anomalous en- 
tropy still present at low temperatures. 

Our proposal for the physics of Er2Ti207 is essentially 
equivalent to previous studies of multiplc-q magnetism 
in face-centre-cubic magnets which contain most of 
the fundamental ideas and concepts. We organise the 
paper into three main sections; firstly a discussion of the 
main experiments and how they lead us to our proposed 
ground-state, secondly a theoretical model that can be 



used to understand the majority of the experiments and 
thirdly a semi-classical spinwave analysis of our model. 
Our model is clearly only approximately solved which 
causes clear 'gaps' in our reasoning. Finally we conclude. 



II. EXPERIMENTS AND INTERPRETATION 

We commence with the specific heat measurements [al , 
which have been repeated as a function of fieldfl, [sl [l6j. 
although the original measurements are all that we re- 
quire. These specific heat measurements lead directly 
to three crucial physical facts. Firstly there is a clear 
phase transition near in energy to the natural dipolar 
energy-scale Secondly there are only two states per 
atom involved in this low temperature phase transition, 
a Kramer's doublet, and so we need only deal with a 
pseudo-spin half. Thirdly the relevant low temperature 
magnetic contribution is gapless and the observed power- 
law is consistent with a magnon dispersion that goes lin- 
early to zero at isolated points in reciprocal-space. It is 
this third point which is totally unexpected: spin-orbit 
interactions provide a coupling between the spin direc- 
tion and the lattice which almost invariably lead to a 
spinwave gap. In transition metals the spin-orbit interac- 
tion is weak in comparison to the exchange and so the 
gap is small, but in rare-earth magnets this interaction is 
usually huge and any associated low energy mode is ab- 
sent. The only likely failure for this argument is when the 
atom has a pure isotropic spin, gadolinium for example, 
where this anisotropy might be expected to be absent. 
Erbium would be expected to have a huge anisotropy 
gap. Sometimes low energy modes are seen (although 
not vanishing), in CeAs for example [isj. however the in- 
terpretation for these modes is different^. These special 
modes occur in multiple-q magnetism, only found in frus- 
trated systems whose magnetism breaks the point-group 
symmetry of the underlying crystal. The modes are as- 
sociated with transferal of spin-density between distinct 
Bragg spots and correspond to local spin rotations which 
are longitudinal antiferromagnetic in character. We are 
proposing that it is one of these special modes that has 
gone soft in Er2Ti207. 

The most important experiments for our development 
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are those involving neutrons. There is a powder diffrac- 
tion experiment [ll , a single-crystal diffraction experi- 
ment in a restricted plane but with critical inelastic 
measurements a fairly complete single-crystal diffrac- 
tion experiment i3|] and a polarised neutron study JJJ . We 
will deal with the elastic scattering first and then consider 
the inelastic separately. Elementary investigation of neu- 
tron scattering data involves three critical elements: the 
structure- factor, the form-factor and the orientational- 
factor. The structure-factor variation stems from the dif- 
ferent moment directions occurring on distinct atoms in 
the magnetic unit cell, and for us is fairly definitive. The 
form- factor arises from the orbital spread of the magnetic 
moment about the site of the nucleus, and is both a com- 
plication and a concern. The magnetic erbium core is in a 
very complicated state, dominated by Hund's rules, it has 
a maximal total angular momentum which is parallel to 
its maximal spin to provide J = 15/2. The sixteen states 
are crystal field split to a relevant pseudo-spin doublet. 
The additional crystal-field doublets are visible using in- 
elastic neutron scattering but are at high energies in 
comparison to the inate energy of the phase transition. 
Indeed, the specific-heat measurements yield almost pre- 
cisely this pseudo-spin degeneracy before the higher en- 
ergy states impinge These complex low energy states 
have a shape that will distort the form- factor [2(?|. We 
will ignore this distortion and presume that the form- 
factor is spherically symmetric. This issue is usually not 
a worry, because domain averaging reinstates the point- 
symmetry, however our proposed magnetism breaks that 
point-symmetry and would highlight these effects. The 
orientational-factor (that stems from the fact that neu- 
trons can only scatter from moments perpendicular to 
the direction of momentum transfer) is very physically 
informative and marks the start of our analysis. 

In simple terms, for any set of spots with the same 
structure-factor, one looks for the smallest (form-factor 
compensated) spot. The magnetic moment associated 
with that structure-factor is then essentially parallel 
to the reciprocal-lattice-vector, G, of that spot. The 
neutron-scattering data then provides the following in- 
formation. The closest spots to the origin, (111), (111), 
(111) (111), are all quite small. The second-closest spots 
to the origin, (200), (020) and (002) are all very small, 
while the third-closest spots, (022), (202) and (220), are 
dominant. The second observation is crucial, since it tells 
us that the three independent components of spin-density 
that characterise the magnetism are parallel to the Carte- 
sian directions. The states are described mathematically 
by 

m(R) =1 m I [e*-^^"sin6'cos0x 

where R is the position of the atom in natural units 
(corresponding to the underlying face-centre-cubic Bra- 
vais lattice) and pictorially by FigH] To limit to the 
depicted magnetism, we also require that there is no 
ferromagnetism in the system, which is controlled by 




FIG. 1: The spin-states consistent with vanishing Bragg spots 
(200), (020) and (002). 



the almost vanishing (222) Bragg spot. Technically, 
these states are equivalent to Type I face-centre-cubic 
antiferromagnetism,4], but interestingly, the pyrochlore 
magnetism is easier to deduce than that of face-centre- 
cubic. The key observation is that the (111) spots pro- 
vide complementary information: this particular spot 
provides a phase which splits the lattice up into planes 
perpendicular to the (111) direction, each of which has 
uniform phase. These planes happen to be alternating 
Kagomc and sparse triangular planes with alternating 
signs. The key observation is that these triangular planes 
make up one of the four sublattices, while the Kagome 
planes make up the other three. If the system is an an- 
tiferromagnet, then the structure factor for this Bragg 
spot should be proportional to the spin direction on the 
triangular sublattice. The experimental observation that 
the (111) Bragg spot is quite small then provides the in- 
formation that the spin on site zero of Fig[l] is almost 
parallel to the (111) direction. This clearly contradicts 
the endemic idea that the spins are perpendicular to the 
natural crystallographic directions[l-3, 5, 14-16I. [l9|. 

The neutron experiments also clearly exhibit distinct 
magnetic domains [l|-|3|. This is clearly inconsistent with 
the idea that the state is a pure triple-q state (spins 
pointing towards the natural crystallographic directions, 
defined by cos 9 = ) , because this state preserves the 
cubic symmetry. The magnetic field experiments provide 
us with an essentially unique state. The experimental ev- 
idence was initially conflicting; the first measurement [l| 
tried to imply domains from a jump in intensity for a 
single Bragg spot when a field was applied. Domains 
correspond to redistribution of Bragg intensity, and so 
if one peak grows then another should reduce, whereas a 
change in a single Bragg spot should signal a fundamental 
change in magnetism and not a change in domain struc- 
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ture. This experimental result was therefore anomalous 
until the single-crystal measurement Q indicated an un- 
expected and unexplained diffuse peak that enlarged the 
Bragg spot with the inclusion of the field. This single- 
crystal work[2] docs show clear indications of domains, 
but because not all of the magnetic Bragg scattering was 
observed it is not unambiguous. The final fairly complete 
measurement [3] measures all three crucial peaks indepen- 
dently and they are seen to separate in an intensity pre- 
serving way once the field is applied. We can fit the three 
intensities to our Eq.(IT]) and we find that 
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although no crystallographic significance should be at- 
tributed to our approximation. This direction is quite 
close to -^(1,1,2), which also holds no special signifi- 
cance, unlike the previously proposed, -^(1, 1, 2), which 
is perpendicular to the natural crystallographic direc- 
tions. We have also analysed the (111) spot [2] which 
appears consistent with this choice, subject to our crude 
choice of form-factor, a pure hydrogenic /-state. The 
state we propose has tetragonal magnetic symmetry and 
lies on the path that smootlily connects the pure triple-q 
state to the single-q state |4j. 

Finally, we consider the inelastic neutron scattering Q 
together with an ESR experiment (l6j . The low-energy 
gapless mode is clearly visible close to the magnetic 
Bragg spots, and a second almost flat mode is also visi- 
ble. Perhaps the most intriguing observation is that when 
the field is applied these magnons disappear, but reap- 
pear close to other point-group related Bragg spots Jj, 
a subtle domain issue that we do not currently try to 
explain. Note that the position of this Goldstone mode 
is exactly equivalent to that in single-q CeAs[18j or in- 
deed in USbflU or U02[2i|, both of which are thought 
to be pure triple-q states. In parallel to this, the ESR 
experiments [16[ see analogous modes but chase them as 
a function of field. The physical surprise here is that it is 
the flat mode that softens and controls the loss of anti- 
ferromagnetism and not the low-energy Goldstone-mode 
as might naively be presumed. This is further evidence 
that the low energy mode is not of the usual type. 



III. MODELLING 

Modelling rare-earth magnetism properly is compli- 
cated. The dominance of Hund's rules makes the orbital 
nature of the states complex and the local crystal-field 
interactions destroy a simple isotropic-spin picture. In 
previous work the case of Gd9Ti907;12j was considered, 
here we anticipate a pure spin with no orbital complica- 
tions and natural models. In this section we develop a 
phenomenological model to try to describe the observed 
behaviour of the erbium compound, and the natural mod- 
els are severely renormalised. 



One major theoretical complication is the crystal-field 
interaction, which projects the original large-J states 
onto a pseudo-spin half representation that controls 
the low temperature phase transition. This projection 
strongly renormalises the interactions and hence corrupts 
intuition. One can provide some understanding using an 
elementary model for this projection. For a half-integer 
J we can easily show that the projection onto the states 
I J, ± J) along the natural crystallographic directions, de- 
fined by. 
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yields the mapping 
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(4) 



where Jq. are the original angular-momentum operators 
and Sq are the pseudo-spin half operators. We can then 
project any bilinear Hamiltonian onto an effective spin- 
half model. For the case of Er2Ti207 the moment along 
the crystallographic axes is thought to be small [sj and 
in addition the ordered moment is indeed observed to 
be smalHl]. To analyse this system it therefore seems 
more natural to project onto the states | J,^\) using 
the following mapping: 



J„ i-^ J 



(5) 



The role of the crystallographic axes is dominant in the 
first case, Eq.Q, and weakened in the second, Eq.®. At 
bilinear order, our projection onto pseudo-spin is quite 
general, if the two coefficients in Eq.® are made param- 
eters, depending only on the idea that the crystal-field is 
dominated by a single natural direction. 

We develop five natural interactions in our modelling. 
Four stem solely from the crystallographic directions and 
bilinear coupling. The first, Hq, has no pseudo-spin ana- 
logue, and the others. Hi, H2 and iJa, provide the domi- 
nant spin interactions. The fifth interaction, H^, couples 
the spin to the lattice and is crucial to our modelling. 
We elect to derive these interactions by looking at the 
natural energies in our system. 

The first energy we consider is the dipolar energy, this 
is not the strongest interaction as both exchange and 
residual crystal-field are expected to dominate, but it is 
very instructive. Working on a single tetrahedron, in 
the presence of dominant antifcrromagnetic Heisenberg 
interactions, the local dipolar interaction can be shown to 
take the form of a linear combination of two contributions 
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and 



This provides a natural model which is minimised by 
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where the coefRcient of Hq is six times that of for 
the local interaction. Using Madelung techniques [Q, it 
can also be shown that the dipole interaction, restricted 
to the q=0 subspace, also takes the same form (for the 
static issues) but with Hq slightly more important. We 
also elect to incorporate the interaction 
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(8) 



which is generated by our projections. Treating these 
Hamiltonians as though they were initially functions of 
Jq, applying the maximal spin projector then yields 



Ha ^ IQJ^m Hn ^ AJ^Hn m ^ AJ^H 



(9) 



All that remains for this case is Hq and H^. Although 
the Hamiltonian Hq is crucial for the original angular mo- 
mentum, once we arrive at our spin-half pseudo-spin, this 
interaction reduces to a constant and can be ignored. As 
we only look at semi-classical solutions to our Hamiltoni- 
ans, we are forced to keep Hq active, although technically 
we need only the single Hamiltonian, i/a, for this limit. 
Although Hq is lost, there is still effectively a residual 
crystal- field interaction, H^ , promoting the natural crys- 
tallographic directions, it is just no longer local. Spin- 
ice is thought to be well represented by the maximal 
projected pseudo-spin. The spin component parallel to 
the crystallographic direction, then a conserved 

quantity and the pseudo-spin always points in these di- 
rections. The residual interaction energy on a tetrahe- 
dron is described by H^ and is minimised by a 'two-in 
and two-out' configuration. 

We will now consider the minimal spin-projector 
model. It is slightly more complicated to analyse but 
provides 



^ ( J{J + 1) + i") ^4 - 4 (j(J + 1)-^]H3 
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The H3 contribution now appears with the opposite sign. 
We observe that in the large J limit Eg. dTUl) becomes 
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H. - AH. = 
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There are many solutions to this, including the natu- 
ral dipolar spirals (as found in Gd2Sn2 07tll]), the cubic 
symmetric triple-q state and appropriate superpositions. 

We now move on to exchange interactions. Suscepti- 
bility experiments suggest that there is a fairly strong 
antiferromagnetic inter action |23|, and for a single tetra- 
hedron this is 
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(13) 



defining a total-spin, T. Although we have chosen to use 
a pure isotropic interaction, the true interaction probably 
depends on the details of the occupied orbitals and is 
consequently likely to be more complicated. Employing 
our maximal spin projection this reduces to 



Hi ^ y J2 {Hq - H3) 



(14) 



and, as expected, conflicts with 'two-in and two-out'. In 
a similar way to the previous interaction we can consider 
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which is the large- J limit. This is minimised by 
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(16) 



providing another natural Hamiltonian. Again there are 
many solutions to this, including all spins being both 
parallel and perpendicular to the natural crystallographic 
directions. 
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Applying the minimal spin projection to Hi provides 



J' 
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where 
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(18) 



Physically we see that the Heisenberg interaction remains 
dominant, but there are also terms equivalent to the max- 
imal spin projector and additionally a new interaction, 
H2-, that would vanish for an antiferromagnet. Note that 
the full dipolar interaction, restricted to a single tetrahe- 
dron, is proportional (up to a constant) to 



Hi + 6Ho - 3H2 



7 



Hi 



(19) 



and so all these interactions are 'expected'. 

Our dipolar arguments naturally lead to Hq, Hi, H2, 
H4 and (through pseudo-spin projection) H3, but our 
use of the Heisenberg interaction is less natural. There 
are likely to be uncontrolled orbitally driven distortions 
to the exchange process in this material. We are con- 
sequently forced into treating our Hamiltonian as phe- 
nomenological. Note that if for each exchange bond we 
simply enhance the bond when the spins are in the Carte- 
sian plane of the bond, at the expense of when the spins 
are perpendicular to the plane, then we generate 



H = {l-S) S^S^ + S^S^ 
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(20) 



as desired in the next section. The spin-orbit interaction 
can easily bias the exchange, which is mediated by orbital 
overlaps, in this way. 

One experimental observation that is problematic for 
the state currently proposed in the literature is the 
severely reduced moment. If the spins ordered in the 
proposed directions, perpendicular to the natural crystal- 
lographic axes, then one would expect to observe essen- 
tially the full moment, whereas a modest ordered moment 
fits the neutron Bragg data. The previously proposed ex- 
planation for this small moment is quantum fiuctuations, 
but our prediction has different physics. We suggest that 



the spins align quite close to the natural crystallographic 
directions, where the pseudo-spin states have a modest 
observable moment, as suggested by Eq.(IS]), in agreement 
with the experiments. 



IV. SPINWAVES 



We elect to use the Hamiltonian 



H = jY,[Hi-SHi~ rjHo + ^{Ho - H^)] 



(21) 



where t labels the tetrahedra, J is the natural energy 
scale, 6, rj and ^ are scale-free parameters. We choose to 
use positive values for these parameters (to stabilise the 
multiple-q states) and consequently one should probably 
view this Hamiltonian as phenomenological. The resid- 
ual crystal-field interaction can account for 77 and ^, but 
the choice of 5 is physically opposite from the presumed 
dipolar source for this interaction. 

This unexpected 'failure' of the dipolar interaction has 
direct experimental consequences. In real-space, point- 
size moments on a lattice have a dipolar energy propor- 
tional to 
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where V{X) is the Coulomb potential between all lattice 
sites but vanishes as X 1— >■ 0. In reciprocal space this 
transforms to 



k G 



(k + G).Sk |2 1/(1 k + G I) (23) 



where V{K) is the Fourier transform of V{X). Now since 
V{X) is best represented as the Coulomb interaction 
multiplied by an increasing function, V{K) corresponds 
to the Coulomb interaction multiplied by a reducing func- 
tion. To minimise the dipolar interaction we should min- 
imise the contributions from the Bragg spots, G, clos- 
est to the origin. This means that the spin-density for 
these spots should be perpendicular to their reciprocal- 
space position. Interestingly, because neutrons scatter 
with essentially an equivalent dipolar interaction, this 
means that these spots should be maximal losing nothing 
from the orientational factor in neutron scattering fl^] . In 
Gd2Sn207, and to a slightly lesser extent in Gd2Ti207, 
these Bragg spots are indeed maximal[llj, whereas in 
Er2Ti207 they are either small or almost vanishing. This 
tells us that the dipolar energy is intrinsically frustrated 
and explains why we are forced into using the apparently 
unphysical sign for H/^. 

Since the system exhibits long-range order, we elect to 
solve our Hamiltonian in the semi-classical limit to assess 
the expected style of order and associated magnons. We 
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employ the Holstein-PrimakofF transforniation[2E 



[|]; (5t + h) 

S-b^ 



(24) 



with the z-axis oriented along the appropriate quantisa- 
tion direction for each spin. We start out solving the case 
77=0=^ to find that the classical solution of Fig [T] denotes 
the ground-state manifold. This then provides 
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S2 = [x y z] 
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for any particular tetrahedron. Bloch transforming pro- 
vides the excitations as 
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where we are using the sublattice degrees of freedom as 
an implicit vector index and where denotes half of 
reciprocal-space, having divided out inversion. The ma- 
trices are explicitly 
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where 

Cq±/3 = cos(fca ± /C/}) (32) 

We can now diagonalise this Hamiltonian with a Bogoli- 
ubov transformation to construct the spinwave spectrum, 
i?kQ • In FigUl we look at the case where 5 is very small. 
There are two almost flat bands on the energy scale of 
2JSy/S. These control the original degeneracy of the 
Heisenberg model which has been lifted by the inclusion 
of H4. In Fig|3]we take the case of a small but relevant 
S. We see that the flat bands have now developed a little 
dispersion and are close to the highest energy excitations. 
This plot has much in common with the experiments Q, 
but also has some clear problems. The low energy mode 
is clearly visible, as is the almost flat high energy mode, 
but instead of a linear dispersion at low energy we find a 
quadratic dispersion. This is because our model is actu- 
ally degenerate across all the multiple-q states of FiglH 
but if the particular ground-state were picked out by 
the energetics then the spinwave would be expected to 



harden. 

We are claiming that the low energy mode observed 
in the experiment is a rare form of longitudinal spinwave 
and so we now take time to explain this idea. We start 
out with an isotopic single-q state, 9=0, oriented along 
the z-axis. The magnetic scattering would then appear 
at (220), but the (022) and (202) spots would be absent. 
The usual Goldstone modes would be expected close to 
(002) and analogous spots on the body-centre-cubic su- 
perlattice generated by {(222), (222), (222)}, but there 
would usually be a spinwave gap at (200) and (020), as 
these points are non-magnetic and unrelated to (002) by 
any magnetic symmetry. The usual Goldstone modes are 
transverse; they are associated with the original state be- 
ing tilted off axis and correspond to a small amount of 
magnetism at the same q-point but in a perpendicular 
direction. The initially gapped modes at (200) and (020) 
correspond to small magnetic distortions of the other two 
styles of point-group related magnetism which are cur- 
rently not present in the system. Physically they amount 
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FIG. 2: The four spinwave branches for the case 5=0.0001 in 
units of 2JS, parallel to one of the Cartesian axes for the case 
corresponding to the experimentally proposed ground-state. 

Spinwave Dispersion 




FIG. 3: The four spinwave branches for the case 5=0.2 in 
units of 2JS, parallel to one of the Cartesian axes for the case 
corresponding to the experimentally proposed ground-state. 



ical idea is then that as the angle d smoothly rotates, 
from zero at the single-q state to cos 8 = in the cu- 
bic triple-q state, the mode remains at zero energy and 

Spinwave Dispersion 



FIG. 4: The four spinwave branches for the case 5=0.2, and 
ry=0.2, which stabalises the cubic triple-q state (<t> = \ and 
Q — Q), in units of 2JS, parallel to one of the Cartesian axes. 



then generates a gap again in the cubic state. Inchiding 
the parameter 77 into the spinwave calculation then sta- 
bilises the cubic triple-q state and provides the gapped 
spin- wave dispersion of Fig|4l 

We now take time to offer an exactly solvable model 
which highlights some of the issues 



(33) 



This model also has the multiple-q states of FigH] as its 
ground-state manifold, but it has the rather elementary 
dynamics controlled by 



H = 2JS 



5+^Z*Z 



■^ZZC\^ 



(34) 



where 



to small changes in the and of FigHJ The softening of 
these modes then corresponds to a phase transition from 
the original single-q state into a multiple-q state. The 
incoming state naturally uses 4> to become a double-q 
magnet or 6 to become a triple-q magnet In a nor- 
mal rare-earth material the spin-orbit interaction opens 
a huge gap in the standard transverse mode, but the lon- 
gitudinal mode, which links to the other magnetic states 
can be found at low energies, as in CeAs^8|]. The phys- 



Ck = 



and 



3 


~Cy + z 




Cx+y 


-Cy + z 


3 


Cx — y 


Cz — x 






3 


— Cy^Z 


— C'x + y 


~Cz-X 


-Cy-z 


3 



(35) 



Z = cos 6 sin (f> + i cos cj) + cos 6 cos (j>~i sin cj) — sin 6 (36) 
The spinwave dispersion is given by 
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ka 



2JS 



(sin 9 sin ( 



cos 9)^ + (cos 6* — sin cos (j)f + (sin 9 cos c 



sm u sm ( 



(37) 



where 



Fluctuation Energy 



ek„ = 4,4,2 ± [1 + 37k] 



(38) 



and 7k is the face-centre-cubic structure factor for the 
underlying periodicity. The quantity eka is essentially 
the pyrochlore lattice structure factor, exhibiting a gap- 
less band and two flat bands describing the degener- 
acy (clear at finite energy). This model (ignoring the 
quadratic dependence at low energy) is enough to de- 
scribe the observed inelastic neutron scattering and 
hence this experiment does not clearly shed much light 
on which multiple-q state might actually be stable. 

Although all our multiple-q states are classi- 
cally degenerate, the quantum fluctuations lift this 
degeneracyd^ and tend to stabilise collinear states. For 
the current exactly solvable model we can calculate this 
quantum fluctuation energy analytically. At zero tem- 
perature we find that the fluctuation energy per spin, Q, 



Q = "iJS 



24 



Z? 



5 S + 
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where 



(39) 



= 3 — (sin 9 cos cj) + sin 9 sin (p + cos ( 
(sin 9 sin (f> — cos 9)'^ + (cos ( 



+ (sin 9 cos (p — sin 9 sin ( 



- sm U cos ( 

\2 



(40) 



This fluctuation energy vanishes for the cubic triple- q 
state and is minimised when the spins are perpendicular 
to the natural crystallographic directions. This exactly 
solvable model has very anisotropic interactions which do 
not promote standard collinear states. If we accept that 
our spin is pseudo-spin half and that the quantum fluctu- 
ation energy is of the same order as the classical energy, 
then we can now attempt a physical explanation for the 
system. If there is a classical energy promoting the cubic 
triple-q state, such as Hq or H3, then this could balance 
the quantum fluctuation energy at the observed exper- 
imental state. However, careful analysis of the current 
model shows that the state jumps discontinuously from 
the cubic state to one with the spins perpendicular to 
the crystallographic directions and so the current model 
is too simplistic. 

Returning to our more physical model. Hi — SH^, we 
can calculate the quantum fluctuation energy as 



Q ^ JS \ -1 -35 



AN 



(41) 




FIG. 5: The quantum fluctuation energy for the case 5=0.2, 
as a function of (in units of |-)with <?f>=-| (solid curve). We 
have included a small rj to achieve the blue curve. 



at zero temperature, and we offer an example in FiglS] 
As expected, the ground-state is now the single-q state 
(black curve). We next add a carefully scaled classical 
contribution favouring the cubic state and now an in- 
termediate state appears relatively stable (blue curve). 
However the energetics are not very 'stable' and rather 
than continuing to evolve smoothly, the ground-state 
tends to make a small discontinuous jump when close 
to the cubic state, to the cubic state. 

Having analysed a particular case, we now move on to 
generalities. Our complete Hamiltonian is 



(42) 



ka 



where is a local interaction which lifts the azimuthal 
symmetry and provides a crystal-field preference appro- 
priate to the weak hexagonal symmetry around the nat- 
ural crystallographic directions. We can semi-classically 
solve such an interaction in general. There are four nat- 
ural ground-states: the dipolar spirals [ll]|, a generalised 
'two-in and two-out' state (typical for spin-icejj), the 
pure triple-q statejlj and finally a multiple-q state with 
the spins perpendicular to the natural crystallographic 
directions, but pointing along one of the twelve natu- 
ral directions allowed by the weak hexagonal interaction. 
The key observation is that in all of these phases there 
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is a natural gap in the spinwave spectrum. Our generic 
Hamiltonian is very anisotropic and that generates the 
expected loss of continuous symmetry and a consequent 
gap. The observed gapless state is truly unexpected. 

There is no easy way to eliminate the gap in the spiral 
and 'two-in and two-out' phases, but since Er2Ti207 has 
a multiple-q ground-state this is irrelevant. One way to 
close the gap is to omit the hexagonal crystal-field inter- 
action, e I— 0. This is the current core of the picture 
prevalent in the literature [ll-lsl, [T5| . If it were exact, then 
one could use a magnetic field to access this degeneracy 
but the current experiments point to the particular az- 
imuthal angle being pinned and nonreactive to a field. A 
weak hexagonal field provides a small gap. However it 
is possible that the observed linear dispersion, and spe- 
cific heat, have not been accessed at low enough ener- 
gies, and temperatures respectively, to observe this gap. 
Another way to close the gap, our explanation, is that 
the multiple-q degeneracy is active. We would leave the 
hexagonal interaction, which stabilises (p = j, but omit 
the excess local crystal-field term, 77 i—)- 0. The inclusion 
of the hexagonal interaction removes one inconsistency 
in our model; the quadratic dispersion at low energy be- 
comes linear as required, and our continuous degeneracy 
is in the angle 6 and not (p. We believe that once the 
particular choice of 9 is chosen by some collective com- 
promise, the gaplessness still remains. Since there is no 
such collective competition in our semi-classical analysis, 
we have resorted to suggesting that quantum fluctuations 
could provide such an asymmetric 9, without offering any 
rigorous calculations to prove this. 

V. CONCLUSIONS 

Rare-earth magnets are usually very anisotropic with 
a gap to magnetic excitations. Er2Ti207 is an unusual 
example of a gapless rare-earth magnet. The experimen- 
tal evidence points towards a multiple-q antiferromag- 
net where the spins do not point along crystallograph- 
ically natural directions. The low energy mode then 
smoothly varies this incommensurate angle. The single- 
domain single-crystal measurements indicate that the an- 
gles (f>—j and (the incommensurate) ^' ~ f are selected. 

The observed spinwave spectrum is well represented by 
the natural pyrochlore structure factor and hence does 
not restrict the model much, except that the spectrum is 
gapless and so a continuous degeneracy is required. The 
observed orientational factors in the neutron scattering 
are contradictory to the expectations of the dipolar in- 
teraction and this leads us to propose the model with the 
opposite sign for in order to flt the experiments. 

Our physical insight arises from the balance between 
quantum fluctuations and classical interactions. Our 
Hamiltonian, H1-5H4, is classically degenerate, with the 
states of FigH] as the ground-state manifold. The quan- 



tum fluctuations stabilise the single-q state, and because 
the pyrochlore lattice is highly frustrated the quantum 
fluctuation energy is expected to be large. Due to the 
strong spin-orbit and crystal-fleld effects there are ad- 
ditional interactions which happen to favour the natural 
crystallographic directions. If we imagine increasing such 
an interaction smoothly from zero, at some strength it 
will overturn the quantum fluctuation energy and drive 
the system into a multiple-q state. We are proposing 
that this transition would be second-order and would 
pass through the intermediate gapless phases where 9 
varies from zero to O, with Er2Ti20i part- way through 
this process at the observed value of 6* ^ ^ . 

Our proposal requires a multiple-q state with an un- 
natural angle 9 which is not clearly derived nor convinc- 
ingly explained. The current picture in the literature has 
the natural angles ^=f and (j>—j and so would be ex- 
pected to have a gap. This gap needs to be too small 
to be observed. This state should also have a dominant 
(111) spot, as this Bragg spot structure factor is per- 
pendicular to its momentum transfer, whereas the spot 
is observed to be small. The only weakness in our ar- 
gument is the form-factor, which would have to have a 
severe and unexpected orientational dependence. 

There are now two plausible states that will flt the 
experiments: the original ideaflS?], with spins perpendic- 
ular to the natural crystallographic directions, and our 
proposal with a non-symmetric angle 9. How can we 
decide between them experimentally? The original pro- 
posal should exhibit a gap, controlled by the hexagonal 
crystal field, which can be sought. If this gap is tiny, 
then an external magnetic field should be able to rotate 
the angle (fi and dramatically change the magnetic Bragg 
scattering, with the ratio of intensities changing from 
(1,1,4) to (3,3,0) for the three main Bragg spots. One 
would also need an explanation for why the (lll)-type 
Bragg spots were so small, cither from dramatic form- 
factor dependence or from some other as yet unresolved 
experimental issue. Our proposal has rather different is- 
sues. Since we fit the known data the experiments are 
not too pertinent, it is the theory that is critical. We 
need a model that stabilises the unsymmetric state and 
we need the state to remain gapless, both issues having 
been side-stepped in this article. 

Our proposed model is a pseudo-spin half model on the 
pyrochlore lattice but our semi-classical magnetic tech- 
niques are not sophisticated enough to predict our ex- 
pectations. Better theoretical analysis of the quantum 
mechanical model is required to make a more detailed 
comparison with experiment. 
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